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Using ★-calculus on the dual of the Borchers-Uhlmann algebra endowed with a combina- 
torial co-product, we develop a method to calculate a unitary transformation relating the 
GNS representations of a non-quasifree and a quasifree state of the free hermitian scalar 
field. The motivation for such an analysis and a further result is the fact that a unitary 
transformation of this kind arises naturally in scattering theory on non-stationary back- 
grounds. Indeed, employing the perturbation theory of the Yang-Feldman equations with 
a free CCR field in a quasifree state as an initial condition and making use of extended 
Feynman graphs, we are able to calculate the Wightman functions of the interacting and 
outgoing fields in a (j>^-iheoTy on arbitrary curved spacetimes. A further examination then 
reveals two major features of the aforementioned theory: firstly, the interacting Wight- 
man functions fulfil the basic axioms of hermiticity, invariance, spectrality (on stationary 
spacetimes), perturbative positivity, and locality. Secondly, the outgoing field is free and 
fulfils the CCR, but is in general not in a quasifree state in the case of a non-stationary 
spacetime. In order to obtain a sensible particle picture for the outgoing field and, hence, 
a description of the scattering process in terms of particles (in asymptotically flat space- 
times), it is thus necessary to compute a unitary transformation of the abovementioned 
type. 
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I. INTRODUCTION 

In the theory of quantum fields on curved spacetimes, it has been realised quite early that in- 
teractions contribute significantly to particle production due to non-stationary gravitational forces 
The scattering theory of interacting fields on physical backgrounds has thus been in 
the focus of attention. The picture that has been drawn is that this scattering behavior is es- 
sentially encoded in two S-matrices, one being a generalisation of the known scattering matrix 
in fiat spacetime that can be calculated by means of reduction formulae similar to the Lehmann- 
Symanzik-Zimmermann (LSZ) method and the other being a Bogoliubov transfor mation between 
"in"- and "out" -representations that are both quasifree but not the same [BiTaSO*]. The former of 
these S-matrices was assumed to describe the scattering of incoming particles into outgoing parti- 
cles via an interaction with both gravity and matter fields, while the latter, and already well-known 
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one 



Wa79f |. was interpreted as describing the particle production by means of the gravitational 



field alone. 

With some exceptions, e.g., BiFo79l | . that, however, have not been very precise on th e topic 



of asyinptotic conditions either, most of the works within the above described framework BirSd . 
BiPaSl iBDFSd . liiT^ have thus assumed asymptotic conditions for the interacting field in the 



Heisenberg picture (j){x) —>■ (/)™/°"*(x) as ^oo with (/)™(x) = a{x) + a\x) and (j)°^^{x) = b{x) + 

b^{x), where a(x), a^{x) and b{x), b^{x) are annihilation and creation fields fulfilling a(x)0™ = 
b{x)Q°^^ = for suitable "in"- and "out"-Fock, i.e., quasifree and pure, vacuum states Jl^"^ and 
Qout respectively. It seems, however, that the consistency of this type of asymptotic condition has 
never been checked. In fact, once the "in"-state for the field has been specified to be a certain, 
say, quasifree state, the field theory should be determined by this initial condition, the canonical 
commutation relations (CCR), and the equations of motion . Prope rties of "out" -fields should 
thus have to be derived and cannot be assumed ab initio. In [GoTa03l | . this supposition has been 
investigated within a toy model for quantum field theory on curved spacetimes which has indeed 
failed to reproduce the assumption of quasifreeness of the outgoing state and has thus underlined 
the relevance of non-quasifree states for free fields. 

In the present work, we show that this result is not an artefact of the mentioned toy model, but 
occurs genuinely in (scalar) quantum field theory on non-stationary spacetimes. In fact, we are able 
to show that the cosmological time scale under which these finding result in non-neglegible particle 
production is roughly one millonth of a second - early enough not to conflict with well established 
physical theories like nucleosynthesis, but late enough to be relevant for physical models of the 
very early universe. 

This raises the question how to physically interpret non-quasifree representations of the CCR, 
i.e., how to compute the particle content of a non-quasifree state. In the case of unitary equivalence 
to a quasifree representation and for real scalar quantum flelds, we provide a complete solution to 
this problem using a ^-product and the associated calculus on the dual of the Borchers-Uhlmann 
algebra endowed with a combinatorial co-product. On a Lorentz manifold that is asymptoti- 
cally flat at early and late times, one can then take the distinguished quasifree "in"- and "out"- 
representations as a reference and calculate the "out" -particle content of the scattered non-quasifree 
state to flnally obtain information on the relative particle production due to interaction with both 
matter flelds and gravitation. 

In order to be able to access the question whether the outgoing fleld is in a quasifree represen- 
tation, one has to explicitely calculate its unsymmetric vacuum expectation values (VEVs), i.e., 
Wightman functions. This is not possible via symmetric scattering amplitudes of LSZ-type or 
Feynman path integral formalisms, one needs a standalone way to perturbatively compute Wight- 
man functions instead. Such a formalism is readily available on Min kowski spacetime, ii amely, 
the sectorised Feynman graph formalism of Ostendorf and Steinmann Ost84l . ISte93l . ISteOCll | . This 
formalism successfully abandons the CCR and asymptotic conditions, but requires the standard 
spectral properties. It is thus not well-suited for our purposes since it can only be extended to 
stationary spacetimes, where spectral properties can be formulated, and can not be employed to 
compute Wightman functions of outgoing flelds, but only those of interacting fields. 

After fixing the setting of our work in section |lll we therefore develop an independent formalism 
to calculate the Wightman functions for quantum fields on generic curved spacetimes using a 
quasifree "in" -representation and the perturbation theory of Yang-Feldman equations, c/., section 
mil In the course of this, we introduce a Feynman graph calculus for Wightman functions where 
the external vertices of the Feynman graph have to be labeled properly and the edges have to be 
decorated with arrows symbolising various kinds of propagators. Once the method to compute the 
Wightman functions has been established, their basic properties are then analysed in the sections 
IIVI and IVl Hermiticity, invariance under orthochronous isomorphisms, perturbative positivity, the 
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relation to asymptotic conditions (l){x) x° ±00 in the Heisenberg picture, and the 

spectrum condition (on stationary spacetimes) do not pose major problems. However, a proof of 
locality by means of our Feynman graphs seems to require a lot of combinatorial effort. We thus 
take a slight detour: after assuring equivalence of the Yang-Feldman equations and the expansion 
of the interacting field by means of retarded products, locality follows easily from the Glaser- 
Lehmann-Zimmermann (GLZ)-relations for retarded products. In this context, the power of the 
GLZ-relation is sourced from the fact that they seem to be the most efficient way to encode the 
commutator combinatorics appearing in the proof of locality. 

Once we have shown that the Wightman functions in our framework possess all sensible prop- 
erties we can expect, we proceed in section IVII by verifying the Klein-Gordon equation and the 
CCR for the "out" -field. The proof of the CCR is again best achieved by having recourse to 
retarded products and it is moreover established that the "out" -field representation is genuinely 
non-quasifree. To assess the physical content of such a non-quasifree state, we develop the afore- 
mentioned method to explicitely calculate a unitary transformation between a non-quasifree state 
and a quasifree one, provided it exists, in section IVIIi In section IVIIH we provide a brief out- 
look. Technical and computational details on ★-calculus and retarded products are given in two 
appendices. 

We point out that all calculations are carried through in unrenormalised perturbation theory. 
This is mostly due to the fact that we have little to contribute to the theory of renormalisation on 
physical backgrounds, w hich has recently been established by means of the Epstein-Glaser method 
and microlocal analysis [BFKQGI . iBrFrOd . iHoWaOll . iHoWaOl EoWa03l | . 



II. SETTING 

In this section, we present some basic notation and the field theory setting we work in. To wit, 
we would like to perturbatively calculate Wightman functions of hermitian scalar quantum fields 
on a globally hyberbolic smooth Lorentzian manifold (M, g) in (/>P-theory. That is, our quantum 
fields, operator valued distributions on a Hilbert space, satisfy the formal^ equation 

{D + + kR)(P = -Xclf-^ 

with coupling constant A, scalar curvature R, mass m and □ = VaV", V being the covariant 
derivative associated with g. 

We s tart by introducing the fundamental functions of the theory. Let Gr {Go) be the unique 



BGP07l | retarded (advanced) fundamental solution of the Klein-Gordon operator (□ -|- m +kR), 
i.e., Gr/a are real valued bidistributions on M satisfying 

(□ + + KR)G^ia{x,y) = 6{x,y) 

and suppa; Gr/a{x, y) C {y X ) being the closed causal forward (backward) cone with base- 

point X. 5{x,y) is the delta-distribution associated with g, i.e., fj^j^j^^dgxdgy5{x,y)f{x)g{y) = 
j^jdgxf{x)g{x) for all compactly supported (complex-valued) test functions f,g G V{M) : = 
Cq°{M,C), where dgX = dx, g = det^i, is the canonical volume form associated with g. We 
note that Gr{x,y) = Ga{y,x) and define the antisymmetric bidistribution 

D{x,y) = Gr{x,y) - Ga{x,y). (1) 



^ Recall that we work in unrenormalised formal perturbation theory throughout this paper. 
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Obviously, D fulfils the Klein-Gordon equation in both arguments and vanishes for x _L y, i.e., for 
spacelike separated x and y. 

To calculate the Wightman functions, we need to specify initial conditions for the field (j){x). 
We achieve this by postulating that for large asymptotic times x^ — > =Foo the interacting field 
(f){x) converges to incoming or outgoing fields where we demand that the "in" -field 

satisfies the free Klein-Gordon equation (□ + + Ki?)(/>™ = 0. For space-times {M,g) that are 
"large" enough to allow the fields to disperse quickly enough to become finally non- interacting, 
the aboy ementioned asymptotic conditions are formulated in terms of the Yang-Feldman equations 
YaFe5nl | 



b^°'(x) 



4>{x)=r/°-\x) + {Gr/aj){x), (2) 

where {Gj-ja j){x) stands for (the formal expression) 

Gr/a{x,j) := / dgyGr/a{x,y)j{y) 



M 



and the current j equals —\(ff~^ in our case. 

It is necessary to specify a representation for (/>™(x) (or, equivalently, an algebraic "in" -state 
for the Borchers-Uhlmann algebra of the free scalar field) "by hand" as, in absence of isometries 
and spectral conditions on gener al curv ed manifolds, the Klein-Gordon equation and the OCR are 



not sufficient to fix it uniquely |Wal93l ] . 



This can be accomplished by first choosing propagators D^{x, y), that is, complex valued bidis- 
tributions on M satisfying the Klein-Gordon equation in both arguments and ImD^ = i^-D, 
D~^{x,y) = D'^{y,x) = D~{y,x), such that D := 21!ieD^ is symmetric and constitutes the 
choice of a state. Here, the bar denotes complex comjugation. To sele ct a pure state, we require 



D{f,f) = lmih(.viM) \D{f,hf/b{h,h)\ for all / G V{M), cf., [Wal93l ]. Particularly, this implies 
that is positive, i.e., D^{f, /) > V / G T){M). We furthermore demand that is invariant 
under any existing isometric diffeomorphisms of {M, g) preserving the time direction, which only 
constrains D, as D automatically fulfils this conditi on. For a discussion of the existence of such 
(and even more general) bi-distributions, cf., Wal93l | . 



Before proceeding to select an incoming state, we need to introduce the notion of truncated 
Wightman functions. These are defined via a cluster expansion as 

{n,cp'^^{xi)---cp'^-{xM= n (Sl.r^^{x,,)---r^i[x,^)^l)^, (3) 

/eP(") {ii,-- ,3iW 

where aj G {in, loc, out} and 'p(") is the collection of all partitions of {1, • • • ,n} into disjoint, non- 
empty subsets {ji, ■ ■ ■ with ji < ■ ■ ■ < ji. A quasifree state is characterised by the property to 
have vanishing truncated Wightman functions for all n 7^ 2. 

Let us anticipate at this point that, via the Yang-Feldman equations, both cj}"'^ and </)°"* are 
formal power series in — A with monomials in 0™ as coefficients and can thus formally be understood 
as operators on the same, i.e., the incoming, Hilbert space, of which is a vacuum state as we 
will specify in the following. This motivates our choice to write all Wightman functions as VEVs 
w.r.t. Vt. 



The standard "Fourier" spectrum condition has been successfully replaced by a microlocal spectrum condition 
[BFK9 6] to advance quantum field theory on curved spacetime in many ways. The microlocal spectrum condition 
does, however, not determine a unique state, but only a class of states [Ver92| . 
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We can now finally specify the state for the incoming field as a quasifree state with two-point 
function D~^, i.e., 

{n,rix)riy)n)^ = D+ix,y), (4) 
{Q,(l)'''{xi)---<p'''{xn)nf = forn^2. (5) 

Taking ([I])-© into account, it follows immediately that ([3]) simplifies considerably if we only 
consider Wightman functions of "in" -fields, namely. 



D^{xj^,Xj^) if n is even, 
if n is odd. 



Here, P'^") is the collection of all partitions of {1, • • • ,n} into disjoint subsets containing two ele- 
ments {ji, ^2} with ji < j2, i.e., is a collection of all possible pairings made out of {1, • • • ,n}. 

We now explain how (jj]) and ([5]) determine a particle picture in the remote past. By our 
assumptions, D is a symmetric bidistribution that fulfils the Klein-Gordon equation in both argu- 
ments. Consequently, the solution part Sf of any test function defined as Sf := D f solves the 
Klein-Gordon equation and clearly {Sf,Sh) := D{f,g) constitutes a well-defined inner product 
on the space of complex solutions of the Klein-Gordon equation with compactly supported initial 
data. Let us indicate the completion of this space w.r.t. ( . , . ) hy TC and note that the imaginary 
part ^-D of defines a (C-bilinear) symplectic form S on the space of complex valued solutions 
via J^{Df,Dg) := D{f,g) that extends continuously to Ti. 

Upon comparison with the symplectic form, the inner product then induces a complex structure 
J via {ip, Jx) '■= ^(V'jX) foi' any two solutions. One straightforwardly obtains J* = — J, = — 1, 
and, thus, J = i{K~^ — K~) with the projector on the eigenspace of J with eigenvalue ±1 In 
the following, we call Ti^ := K^Ti. the positive/negative frequency spaces respectively. We note 
that TC^ = Ti^ since Jtp = —itp for tp £ TL'^ . 

Let now be the symmetric Fock space over Ti'^ with Fock- vacuum 0. By a^('i/'), 
ip G 7^"*", we denote the usual creation and annihilation operators on T. We use the convention 
a(V')* = a^Tp)* , ip G W"'", in order to obtain a C-linear definition for a(x), X £ ■ Here, * stands 
for taking the adjoint (neglecting domain questions). Let S± := K^S be the operator that maps 
test functions to the positive/negative frequency solution part. The incoming field can now finally 
be defined as the C-linear operator valued distribution 

r(/) = a(5_/) + at(5+/). (7) 

Furthermore, by |Wal93l . Lemma 3.2.1], (Z)/, Dh) = D{f, h) = T.(Df, Dh) = {Df, JDh) for all test 



functions f,h, from which —JDf = Df fol lows. A pplication to ^ gives = ia{K Df) 

ia' (K^Df) which is the definition given in [ Wal93l ]. 

It follows from the Fock construction, and the properties of D~^ that 

[cP^^ix),riy)]=iDix,y) (8) 

which constitutes that the incoming field fulfils the OCR and closes the specification and analysis 
of the properties of the "in" -field. 

Fixing both (p^^{x) and would over-determine the system, we therefore only employ the 

Yang-Feldman equations ([2]) as a definition of (p°^^ (x) without specifying any further properties of 
it. From ([1]) and ([2]) it follows immediately that 

cp--\x)=cP"^{x) + {Dj){x). (9) 



Non-quasifree states on curved spacetimes 



6 



Furthermore, because both and D fulfil the Klein-Gordon equation, does as well. 

However, we still need to check if the outgoing field fulfils the CCR and determine whether it is in 
a quasifree state or not. 



III. CALCULATION OF THE WIGHTMAN FUNCTIONS 



To evaluate Wightman functions, we will make use of generalised Feynman graphs. In the 
following figures, we draw all graphs in (/)^-theory for simplicity. For the actual calculations, the 
degree p of the (?!)^-theory is irrelevant. We begin developing the graphical calculus by introducing 
the symbols for the propagators of our theory. 

£)+(x,y) = X . < . y = D'(y,x) 

Grix,y) = X . « • y = Ga(y,x) 

D{x,y) = X -^a— y = -D{y,x) 

D{x,y) = X .O* y = D(y,x) 



FIG. 1 Propagators 



D~^{x,y) is being represented by a line with an open arrow, D{x,y) by a line with a closed 
arrow and Gr{x,y) by a line with a double open arrow, the arrows pointing to x. Furthermore, 
D(x,y) is drawn with two arrows pointing apart as shown in figure [TJ 

Next, we introduce a tree expansion for the fields according to the Parisi-Wu method PaWuSl]. 
We expand the fields in powers of the negative coupling constant —A: 



\x) = f2i->^rrAx). (10) 



a=0 

Clearly, 



^''^^^ 1 otherwise. ^ ^ 



We calculate (p'^{x) for a = loc/out recursively using ([2]) and 



p-1 



A" 



E n 



iloc 



I V 



>-l= 



if CT = and 



(x) otherwise. 



(12) 



where A'°= := Gr, A°"* := D. Following the recursion in (fT2]) . we define tree graphs corresponding 
to the summands in p2|) by an induction over a. To fix the initial step, we draw (p^ix), i.e., an "in"- 
field, as a leaf attached to a root corresponding to an external x-vertex. A tree corresponding to 



is drawn by taking p— 1 trees corresponding to perturbative local fields of order ai, 



1 Cp-i 



s.t. ^ (7/ = (T — 1, assembling their roots yi, - ■ ■ , yp-i to form a single internal y- vertex and adding 
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a trunk, i.e., a new line from y to x corresponding to A°(x,2/). Therefore, a tree correspoding 
to (p^ix) has a root corresponding to an external x-vertex, a trunk corresponding to A", several 
branches corresponding to GrS, several leaves corresponding to "in"-fields and a branching points 
corresponding to internal vertices with a total number of p — 1 branches and leaves emerging from 
them. We note that the causal flow (the direction of the G^-arrows) always points to the root. 

We label the different tree components inductively, accounting for the fact that the indices ai 
in (I12p are distinguishable. The initial step of the labelling induction is fixed by defining the label 
of a trunk to be the index of the external vertex variable attached to it. To assign a label to a 
branch/leaf, one takes the label of the branch or trunk the considered branch/leaf emerges from as 
a basis. One then extends it by a dot followed by a number reflecting the position of the considered 
branch/leaf (field) at the actual branching point (in the corresponding current), i.e., the index i 
of (1)^°^ in (|12p. Some examples of trees are displayed in figure [2] for the convenience of the reader. 




FIG. 2 The only possible tree for (f>Q{x), the only possible tree for (/)'°'^(x), and one possible tree for 03"*(a;) 

We shall proceed to consider the Wightman functions of our theory. To compute them, it is 
sufficient to consider only their connected parts, i.e., the truncated Wightman functions. 

In order to calculate the truncated Wightman functions, we will first expand them in powers of 
the coupling constant: 

oo 

{n, (xi) • • • 4^'^- {xnW = Y.^->^r{^, (xi) • • • {xnW,. (is) 

(T=0 

Inserting (jlO|) into the left side of (jl3p and comparing terms of equal order in —A, we get: 

oo 

{n, 4^'^^ (xi) • • • cf- {xn)n)l = Yl "^-i (^i) • • • Mnf. (m) 

fir" : 0-71=0, ^iJi=a 

We know from the tree expansion that, for a fixed a, every 4>'^{x) can be expressed in 
terms of incoming fields convoluted with fundamental functions. Therefore, it follows that 
(ri, i;^"! (xi) • • • (/>"^(x„)r2)-^ can be expressed in terms of Wightman functions of "in"-fields inte- 
grated with additional propagators. Finally, combining ([TT]) . (fT2|) . ([6]) and p^ . we can express 
truncated Wightman functions of arbitrary fields merely in terms of fundamental functions. 

Let us now introduce Feynman graphs corresponding to perturbative n-point Wightman func- 
tions. A Feynman graph of order a consists of n external vertices corresponding to the arguments 
xi, • • • ,Xn and type-indices ai, • • • , of a Wightman function and a internal vertices correspond- 
ing to arbitrary variables in M. The vertices are connected to the remainder of the graph by q lines, 
with q = 1 {q = p) for external (internal) vertices. A line is called an external line if it is connected 
to an external vertex, an internal line otherwise. While Wightman functions correspond to Feyn- 
man graphs, one can show that truncated Wightman functions correspond to connected Feynman 
graphs. We call a Feynman graph with arrows and labels on all lines an extended Feynman graph. 
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On the level of graphs, the resolving of (jl4p via ([6]) corresponds to gluing the leaves of n 
trees with a total order of a together to yield an extended Feynman graph of order a with n 
external vertices. As we calculate truncated Wightman functions, only gluing possibilities that 
yield connected Feynman graphs are allowed. 

To finish describing the gluing process, we need to analyse the gluing lines in more detail. A 
line originating from gluing a pair of two leaves together, i.e., a D^-line, is defined to be labelled 
by combining the leaves' labels to a pair. Starting from the beginning, we compare the two labels 
slot-by-slot until we find a pair of numbers that does not match. The arrow on the D^-line then 
points to the leaf corresponding to the lower of these numbers. 

We have now described how to expand fields to trees that are subsequently assembled to ex- 
tended Feynman graphs (see figure [3] for two examples). 



X2 




FIG. 3 Two possibilities to glue the same set of trees to graphs corresponding to 

but this also works the other way round. To calculate (17, c/)'*^ (xi) • • • (/>""(x,i)r2)^, we draw all 
topologically possible connected Feynman graphs of order a with n fixed external vertices xi, • • • ,Xn 
of type ffli,--- ,a„. We then consider all possibilities to partition each Feynman graph into n 
connected and loop-less subgraphs, i.e., trees and several remaining lines. Each such subgraph 
contains exactly one root Xi of type Oj. A partition is fixed by marking a certain number of lines 
such that the Feynman graph with these lines removed consists of n disconnected tree graphs 
without leaves, arrows and labels and each internal vertex is part of a such a tree (see figured] for 
two examples, where the marked lines are displayed as dotted lines). In this process, the external 
lines connected to external vertices of type a = in always have to be marked. Next, we assign 
labels to all external lines according to the indices of their external vertex variables. Using these 
labels as an initial step and starting from the roots, we "walk up" the trees on the unmarked 
lines and inductively assign labels to all lines emerging from the vertices we pass. As a result, the 
marked lines have two labels, one from each vertex they connect. The inductive dependence of the 
labels on the preceding labels is fixed by the labelling algorithm defined above in the discussion 
of the tree expansion. For each of the possible choices of labels we draw arrows on all lines. The 
type of the arrows on the unmarked lines is chosen according to ai, • • • , a„ and the tree expansion. 
Furthermore, each marked line becomes a D^-line. The direction of the arrow on such a line is 
determined by comparison of the two labels of the line in the manner described in the preceding 
paragraph. 
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FIG. 4 Two possibilities to extend a Feynman graph corresponding to (fJ, (/)'°'^(xi)(/)'°'^(x2)0°"*(a;3)ri)jf 



To obtain the analytical expression corresponding to an extended Feynman graph, we assign 
variables to all internal vertices, write down the propagators corresponding to all lines and then inte- 
grate over all internal vertices. Once we have the analytical expressions, summing over all topologi- 
cally possible Feynman graphs and all possibilities to extend them yields (0, (xi) • • • </>'^" (x„)0) J 
(see figured] for two examples). 

For calc ulations , it is often convenient to drop the labels and replace them by combinatorial 
factors, see |Hac07l ] for a detailed discussion of these issues. 



IV. PROPERTIES OF THE WIGHTMAN FUNCTIONS: INVARIANCE, HERMITICITY, SPECTRAL 
PROPERTY, POSITIVITY, AND THE ASYMPTOTIC CONDITION 

With the means of computing the Wightman functions of our theory at hand, we can continue 
by discussing their fundamental properties in this section. 



In variance under orthochronous isometric diffeomorphisms As we have shown in section Hill the Wight- 
man functions of our theory can be expressed in terms of integrals of products of fundamental func- 
tions. Since we know from section [II] that all fundamental functions are invariant under isometric 
diffeomorphisms preserving the time direction and the integrals contain the canonical volume form 
which is invariant under all isometric diffeomorphisms by definition, invariance of the Wightman 
functions follows immediately. 
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Hermiticity The Wightman functions fulfil Hermiticity if 

(f^,</.»i(xi) = {n,rH^n)---r'{xi)^)- (15) 

Since we can express the Wightman functions in terms of Wightman functions of "in" -fields con- 
voluted with the real valued fundamental functions G^/a and D and since the order of fields in the 
latter Wightman function corresponds to the order of fields in the original Wightman function, it 
is sufficient to prove (fT5]) for oi, • • • , = in. 

We know from ([6]) how to express Wightman functions of incoming fields in terms of D^. The 
complex conjugation of exchanges all (^Xj^, xj,^ ) with D^( Xj2,Xj-^) which obviously corre- 
sponds to reversing the total order of fields in the Wightman function of "in" -fields, thus (|15p 
holds for ai, • • • , a„ = in. 



Spectral condition In general spacetimes, there is no well-defined Fourier transformation, therefore, 
standard spectral conditions can not be formulated. However, in stationary spacetimes, the time 
translations form a one parameter group of isometries and we assume that Fourier transformations 
w.r.t. the time parameter are possible. A well-defined standard spectral condition can thus be 
formulated in that case: one restricts D by requiring that the Fourier transform in the time ar- 
guments defined by the global timelike killing field of n-point Wightman functions of "in" -fields 
vanish if the sums Sj := all positive. Here, Ei is the variable conjugated to the 

l-ih time argument in the VEV. We note that all fundamental functions are invariant under time 
translations by our assumptions, hence, they only depend on time differences of both arguments. 
Therefore, the unitary time translation operator that is obtained from the time translation invari- 
ance of Wightman functions of interacting and outgoing fields via the GNS construction coincides 
with the time translation operator for the "in" -fields, which has positive spectrum by construction. 



Perturbative positivity If we expand Wightman functions perturbatively up to a given order N, we 
can add further terms of order 0{X^^^) to obtain a VEV of fields (p"''^{x) = '}2a=o(~ '^)'^ ^a-i^) 
that act as operator-valued distributions on the incoming Fock space. The VEVs obviously fulfil 
positivity. Thus, the Wightman functions expanded in —A up to an arbitrary but fixed order N 
fulfil positivity up to a 0(A'^"'"^)-term. 



Asymptotic condition In this work, we have used asymptotic conditions given by the Yang-Feldman 
equations ([2]). It is a natural question to ask up to what extent these asymptotic conditions lead 
to the asymptotic conditions in the Heisenberg picture <p{x) as x^ — > ^oo in a given 

foliation Af ~ R x C of M, where C is a Cauchy surface. In a weak sense, the Heisenberg asymptotic 
condition is 



lim 

O^ioo 



{n, i^i)--- Hxj) ■ ■ ■ <A"" {xn)nf - {n, (si) • • • <A°"*/'"(a:,) • • • <A"" {xuM 



0. 



(16) 

Let us consider the case x*^ — > — oo first: In the extended Feynman graph expansion of the left 
hand side, all connected graphs where the tree with root Xj is of order aj = cancel and all 
other graphs survive. Likewise, if x^ — > +oo, we obtain in the expansion into connected extended 
Feynman graphs two times all graphs where the order of the xj tree is larger than zero - once with 
the trunk of the j-th tree being a retarded propagator for the local field and once another graph 
with opposite sign where the trunk is evaluated with D. Using D = Gr — Ga, we see that, in the 
case x^ — > +00, one obtains all extended graphs where the tree with root Xj is of order larger than 



Non-quasifree states on curved spacetimes 



11 



zero and its trunk is evaluated with an advanced propagator, c/., figure [5l We note that in the 
hmit — > ±00, the integral over the vertex variable u is restricted to u in the causal future/past 
of Xj and hence the domain of integration becomes smaller and smaller. An actual proof of the 
vanishing of the left hand side of ()16p requires technical assumptions on the manifold M and on 
the propagators - and hence the state - and we do not want to go into the details now. It 
however seems that these conditions are not much stronger than what is needed to assure that the 
integrals over the vertices in the Feynman graphs exist. 




> 



FIG. 5 Heuristic check of the asymptotic conditions on the level of graphs 



V. PROPERTIES OF THE WIGHTMAN FUNCTIONS: LOCALITY 

To prove locality of the truncated Wightman functions, we have to show that 

{n, (xi) ... (xi), r-^' ■■■r'^ {xnW (17) 

vanishes for Xi _L Xj+i, Vz E {1, ... ,n — 1}, aj = loc if j G {hi + 1}, aj = in/loc/out otherwise. 
For proving this, it is sufficient to show that the interacting field itself is local. Since we are given 
the interacting field as a formal power series in —A, locality of (/> = (/>^°^ has to be proven to each 
order in —A separately. To order a we have 



[m,m]„= E [4>{xUA{y)<r,]. (18) 

CF\+(J2=CF 

For a fixed order, it is hence possible to replace the fields (/'(x)o-i , 4>{x)a2 by their tree expansions 
and compute [(t>{x)-,4>{y)\^ as sums of commutators of single trees that are in effect commutators 
of products of free fields integrated with retarded propagators. Employing Leibniz' rule for the 
commutator results in glueing together two leaves, one from each tree, with a D-propagator. 
One can then hope to obtain an expression which vanishes for spacelike-separated x and y. The 
procedure for o" = 1 in (j)^ -theory is depicted in figure[6l where the last step follows from D = Gr—Ga 
and "telescope cancellations". 

It turns out that one is left with a "Gr--chain" and a "Ga-chain", i.e., a product of re- 
tarded/advanced propagators such that the right slot of one propagator corresponds to the left 
slot of another, namely. 



[0(x),0(?/)]i = 2i J dgXi {Gr{x,xi)Gr{xi,y) - Ga{x,xi)Ga{xi,y)} (p'^'ixi). (19) 

M 

Owing to the causal support properties of the retarded and advanced propagators, we see that 
either x ^ xi ^ y or x ^ xi ^ y , where x ^ y (x ^ y) depicts that x £ Vy . As a result, both the 
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[(j)>-(x),(j)i<"=(y)]j = 



I 





t 



= 2 




= 2i 



1 



+ 2 




+ 2i 



1 



— < I » ■ ■ « i <1— 

X y X y 



= 2i 



J. 



2i 



1 



« I « ■ ■ » I » » 

X y X y 



FIG. 6 The commutator of two interacting fields to first order 



G^-chain and the Ga-chain vanish for x _L y. To generaUse this observation to arbitrary order, we 
prove the equivale nce of t he tree expansion to the expansion into retarded products. We then use 
the GLZ relation |GLZ57l | to prove locality. See also |Hac07t | for a proof up to second loop order 
that works on the level of Feynman graphs and proves locality graph by graph. 

The retarded product -Ri_„(i?o(a;o) I -Bi(xi), . . . , of n + 1 operators Bo{xo), . . . , -B„(x„) 

that are mutually local, i.e., [Bi{xi), Bj{xj)] vanishes for all i, j if Xi _L xj, is defined as 



^1,0 (^o(a:o)) :=^o(a;o), 

Rl,n {Bo{xo) I Bi{xi), . . .,Bn{Xn)) := 

■= (-1)" X] [BTr{n){x.^l^n)), [Bn(n-l){Xn{n-l)), ' ' ' [^ttCI) (a^7r(l)), ^ol^^o)] •••]] la;od2^^(i) :<-r<a^V(n) ' 
vrG5n 

(20) 

where Sn is the permutation group of order n and 1a denotes the characteristic function of the 
set A. Note that renormalisation is necessary for a proper definition at coinciding points. With 
the above definition, Ri^n{Bo{xo) \ Bi{xi), . . . ,Bn{xn)) is both symmetric in the last n slots and 
manifestly vanishing if any of the spacetime positions of the last n operators is not in the causal 
past of xq. From ([20]) it follows straighforwardly that retarded products of a certain order may be 
expressed in terms of retarded products of one order less 

Rl^n {Bo{xq) I -Bi(xi), . . . , Bn{Xn)) = 

n 

= - X] [Bj{xj),Rl;n-l{BQ{xQ) I Bi{xi),.. . . . . ,S„,(x„))] Ix.^^x, Vi6{0,...,n}- (21) 

i=i 

Employin g the ab ove recursion and the Jacobi identity for the commutator, one can prove the GLZ 



relations jGLZ57l | 



i?l,„ {A I C, Sl(xi), . . . , Bn[Xn)) - Rl,n {C \ A, Bi{xi), 5„(x„)) 
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E 

ICN :={!,. ..,n} 



JjBi(Xi) 1 ,Rl,\N\I\ 



C 



jeN\i 



(22) 



The most important feature of retarded products in our context is the fact that the interacting 
field can be expanded in terms of retarded products, where a retarded product of order a encodes 
the complete (— A)'^-contribution of the interacting field, namely, 



E 

(T=0 



(zA)'^ 



a'. 



(23) 



i=l 



with £int(x) := (/>™(x)^/p in our case. A combinatorial proof of (j23p that makes the relation with 
the tree expansion explicit will be given at the end of this section. 
To prove locality in terms of retarded products, we start from 



i=l (T1+(T2=0" 



Cri!(72! 



X [i?i_CTj(0™(x) I £int(xi), . . . , Cint{Xai)), Rl,a2{<P'^{y) I Ant (Xcn+l) , • • . , Cintix^))] 

which, due to the GLZ relations, simplifies to 



where -Ri,n+i {4^^'^{x) \ (/)™(y), £j®") stands in shorthand for 

j WdgXiRi^n+i{4>'^{x) I (y),£int(xi), . . . ,£int(a;n))- 

The result is, as we could have expected from our finger exercise in figure [6l a retarded piece, pre- 
sumably corresponding to a sum of G^-chains, minus an advanced piece, presumably corresponding 
to a sum of Ga-chains. In fact, we will prove this correspondence in appendix [Bl since it will be 
necessary for the proof of the "out"-field OCR. If x _L y, both the retarded and the advanced piece 
vanish due to the causal support properties of the retarded products. 

Let us now proceed to see why the perturbative expansion of the interacting field in terms of 
retarded products is equivalent to the perturbative expansion of the interacting field due to the 
Yang-Feldman equation in terms of tree graphs. Since in case of the expansion in trees the (— A)'^- 
contribution to the interacting field consists of the sum of all possible tree graphs with a branching 
points, we have to show that the (— A)'^ term of the expansion in retarded products corresponds to 
exactly such sum. To achieve this, let us recall how one can inductively obtain all possible trees 
of order a: one starts with the "in"-field and then replaces an "in"-field/leaf by a first order-tree 
a times and in all possible ways, taking care to discard trees occuring multiply. Starting from the 
recursion (j2ip , we shall proceed to understand how it reproduces the aforementioned combinatorial 
procedure. 

To this effect, let us first notice that ()2ip simplifies considerably in the current context, since 
all operators . . . , Bn are of the same type and the labelling of the Xi does not matter due to 
them being integration variables. Hence, the symmmetrisation can be replaced by a factor and the 
(— A)'^-contribution to the interacting field in terms of retarded products reads 
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Y\_dgXi Rl^a{(t>'''{x) I £int(a;i), • • • , Ant (a^a)) 



i=l 



{a-iy. 



ITdgXj [£int(Xa),i?l,<7-l((A"(a;) |£int(a:i),--- ,Ant(a;a-l))] 1 ^rr^^r 

ViS{0,...,<T-l} 



Since all operators appearing in i(<?5'™(a^) I >Cint(a;i), . . . , £int(a;o— i)) are monomials in 

the "in" -field, we know due to the Leibniz rule for the commutator that this retarded 
product is given by sums of products of "in" -fields. As Cintixa) is also a monomial in 
the incoming field and the commutator of two "in" -fields is a C-number, we can compute 
[Cint{x„),Ri^a-i{(p'''{x) I £int(a::i), • • ■ , Cintixa-i))] by means of 



i=l 



m j-l "1 
j=0 1=1 i=j+l 



which holds under the assumption that commutes with all other occuring operators for all 

i and where we have formally written dB^/dB in shorthand for nB'^~^. This formula implies that 
[Cir,t{xa),Ri,a-i{(t>''^{x) \ £int(xi), . . ■ , CirA{xa-i))]'^xa>xiii&{o,...,a-i} Can be computcd by summing 
over all possibilities to replace one incoming field (j)™{xj) in Ri^a-i{(j)^'^'^{x) \ £int(a;i), . . . , £int(a;o-_i)) 
by 



iGr{Xj,Xfj)4>^'^''{XcrY^ 



To account for the Leibniz rule of this procedure, we denote it by means of a formal derivative 
operator (c/., DuFr02l | . where the framework is such that it is a well-defined functional derivative), 



VIZ., 



V{x)r{y) ■■=Gr{-,x)^M^^nir^riy) ■.= Gr{y,x)^^'^'^''^ 



Employing this notation, we can write the interacting field to order a as 



(l){x)a 



i-ir 



[a 

i-i) 



TTT / IT '^3^* [Cint{Xa),Rl.a-l{(l)''^{x)\Cint{xi),...,Cint{Xa-l))]l ^ah^i 
Ij! J Vi6{0,...,a-1} 



cr-l 



(a-1)! 



/ n^^<?2;j P(x^)iii,cr-i((/>"'(a;) I Ant(a^i), • • • , Ant(a:cr-i))l 



i=l 



Vi6{0,...,(T-l} 



where it is understood that the operator expression replacing the "in"-field on which T>{Xf^) cur- 
rently acts has to be inserted exactly at the position where that incoming field has been. We can 
now iterate the recursion ()2ip to eventually obtain 



(p{x) 



-- j WdgXiViXcr) ■ ■ ■'D{xi)<j)''^{x)lx„>->xiyx- 



(24) 
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This puts us in the position to prove formal equivalence of this expression to the one obtained 
from the Yang-Feldman equation in terms of tree graphs. Obviously, both expression are equal in 
zeroth order. Let us assume that they are equal to order a and show how equivalence to order a + 1 
follows inductively. By our assumptions, the integrand of (j24p can be rewritten in a form devoid 
of explicit restrictions on the integration domain for xi, . . . , x^, namely, the sum of the integrands 
of all trees of order a, viz., 

(j){x)a=: j Y\_dgXi laixi, . . . ,Xa). 

With this notation, we have 

ViG{0,...,cT} 

T>{xa+i) acts on by replacing an "in" -field by the integrand of a first order tree graph in all 
tree integrands consists of and in all possible ways. As a result, we obtain integrands of trees 
of order o" + 1, still with no mutual restrictions on the integration variables xi, . . . but with 
the constraint that Xo-+i is causally earlier than all of these spacetimes points. It is clear that 
{xaJ^i)Ia{xi, . . . ,Xa-)lxa+i>iXi\ii(^{o,...,a} coutaius all possiblc tree integrands of order a + 1, the 
question is whether on can rewrite this expression in way that it contains every possible tree 
integrand of order o" + 1 with unit weight and without any restrictions on the integration domain. 

To answer this question, let us divide the tree integrands we straightforwadly obtain by the 
action oi V{xa+i) on . . . , Xq-) into equivalence classes, where two integrands are taken to 

be equivalent if they can be matched by permuting vertex variables. Let us choose an arbitrary 
but fixed equivalence class E and assume that it has m members. Their number implies that 
there are m different possibilities to build the tree graph corresponding to E out of a tree graph 
of order a by replacing a leaf with a first order tree. Hence, the tree graph corresponding to E 
must have m "virgin" vertices, where we call a vertex virgin if it has the maximal number of 
p — 1 "in" -fields attached to it; the m members of E then correspond to the m possible ways 
to remove one of these virgin vertices to obtain a tree of lower order. Now, let us permute the 
vertex variables of all members of E such that they are all equal but have different restrictions 
on the integration domain of xi, . . . , Xo-+i and let V := {xjj , . . . , a^i^} be the resulting integration 
variables of the virgin vertices. Since the virgin vertices are connected to the remainder of the 
tree by retarded propagators, we can discard the integration constrains on the m integrands under 
consideration which are automatically fulfilled due to the causal support properties of Gr- The 
remaining restrictions on the integration domain can only be of the form V 9 Xj ^ Xj Vxj G V. In 
fact, due to the Leibniz rule for I?(xo-+i), all m possible integration constraints of this kind must 
appear. Hence, summing up the m integrands with matching variables but different restrictions 
on the integration domain, we obtain once the same integrand, but without any integration con- 
straints. Since the above described procedure is valid for all equivalence classes of tree integrands, 
/m<^+i nr=jL^ dgXiD{x^j^i) ■ ■ ■ P(xi)(/)™(x)lx^+i^...^xi^a; Corresponds to the sum of all possible trees 
of order a + 1 weighted with unit multiplicity, and the formal equivalence of the expansion of 
the interacting field by means of the retarded products on the one hand and the Yang-Feldman 
equation via Parisi-Wu tree graphs on the other hand is established. 

VI. PROPERTIES OF THE OUTGOING FIELDS 

We now examine the properties of the outgoing fields. 



(j){x)a 



+ 1 



„ cr+1 

/ Y\.dgXi'D{Xa+l)Ia{xi, . . . 
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Klein-Gordon Equation We have already seen in section [IT] that the "out" -fields of the theory fulfil 
the Klein-Gordon equation as is manifest from their definition via the Yang-Feldman equation ([9]). 

Canonical commutation relations To prove the OCR, we have to examine the commutator of two 
outgoing fields to each order in —A separately. Due to the Yang-Feldman equations, we have 

[</.-*(x),0-*(y)]^ = [0(x) - i-X)Gaix,<tf-'),Hy) - i-X)Ga{y,r-')]„ 

= [cl,{x)A{y)L-[m,Ga{y,(^-')]^_,-[Ga{x,clf-^),cl,{y)]^_^+[Ga{x,r^^^ 



iia Ilia m 

(25) 

In zeroth order, only the first term of (j25p contributes and the result is 

[cp°-'{x),r-\y)]^ = [^'\x),(l,'yy)] = iD{x,y). 

To prove OCR for the "out"-field, we thus need to show that (j)°'^^{y)^ vanishes identically 

for fj > 0. In the case a = 1, only the first three terms of (|25p contribute and we have 

h = -i j dgxi {i?i,2 I 0^°(y),£int(xi)) - i?i,2 | ^^'^(x), £int(xi)) } 

M 

= {p-l)i J dgXi {Grix,xi)Grixi,y) - Gaix,xi)Gaixi,y)} (p'^'ixif''^ 

M 

Ih = [(t^%x),Ga{y,{<t^n''-^)] = (p- 1)2 j dgXiD{x,xi)Gr{xi,yW%xif-^ 

M 

[Ga{x,{cl}'r~^),r{y)\={p-l)t j dgXiGa{x,xi)D{xi,yW\xiy-\ 



nil 



M 



Due to "telescope cancellations" by means oi D = Gr — Ga, Ii — Hi — HIi = 0. 

For o" > 1, the structure of cancellations is in principle the same as above, the only difference 
being the possible appearance of propagators "lying inbetween" G^/aix, xi) and Gr/a{xi,y). Such 
terms survive in I„ — II^j — Ilia and have to be cancelled by IV^. To treat such terms, we need 
the following identities, proven in appendix B: 



— j"!" 



M 



(_ ■^o■2 

dgXiGr{x,xi) Y.^{x^y'^r-^Ri,..+i{<t^'\xi)\<t^''{y),i^f2) Hxifaf-' (26) 



da^i E Y.^{xi)i^-^Ri,a,+i{r{x)\r{x,)Xf:,^) c^ixiZf-^ Grixuy), (27) 

^ J2^i=n-lj=0 
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where (p{x)a stands in shorthand for ^ Y\ 4'{x)iJi- Iterating these identies yields that 

^ n\ -^1,^+1 ('^™(^) I '^™(^)' A^") '^^'^ expressed completely in terms of Gr-chains connecting 
X with y. 

Employing the above listed identities, we have for o" > 1 



-if 



p-2 



dgXl ^ ^Ct>{xi) 



"1 C72! 



Ga{x,x,)R,,,,+, {r{y) I r(xi),£f,^^)} c/>(:ei 



\P-2-i 



M 

p-2 



+ / dgXidgX2 ^ ^ </>(xi)*^^0(a;2) 

;V/2 J]o'i='^-2*J=0 



(-0 



- G,(x,Xi)iii,.3+i (r(x2) I Ga{x2,y)] ct>{x2r^f'^ct>{xi)Pf-\ 

where the first summand of the last line corresponds to the case (T2 = in the second line. 

For IIu and Ilia, we again need (j26p . (j27p and, fm'thermore, the general commutator identity 



i=i j=i 



n— Im— Ifc— 1 i— 1 m n 

EEYi^^Yi^Mk^Bi] n n 

k=l 1=1 i=l j=l j=l+l i=k+l 

n—lm—ll—l k—1 n m 

EEn^^n^^i^'^'^'] n n 

k=l 1=1 j=l i=l i=k+l j=l+l 



(28) 



to obtain 
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II„ = j dgXi [(l){x),(l){xiY ^]^_^ 

M 

r "'^ 

M Z;o"i=o— 1 *=o 



M 



J2ai=a-1 i=0 



i-i) 



M 

p-2 



(J3! 



y dgXldgX2 ^ ^ (?:)(X1)^^0(X2) 

X {G,(x,Xi)i?i,.3+i(r(xi)|</>i°(x2),£fnr) - 

-G,(x,xi)i?i,.3+i(r(x2)|r(xi),£f„^^)}0(x2)?72-i^(^^)P-^^^ 



= y" dgXiGa{x,xi)[(t>{xiY-\(t>{y)]^_^ 



M 

p-2 

d,xiG,(x, xi) ^ ci>{x,)i^ [ci>{xii mu <t>i^ir.f~' 
dgxiGa{x,xi) E'^(^i)-i^^{^i.-2+i('A''^(xi)ir(y),/:fr) 

(p-i)./<i,..G„(...,){a(:.,..)-G.(..,.))*(.,)r-; + 

M 

+ dgXidgX2Gaix,xi) 22 z2 ^^^^y<Ji(t>{x2)i2 ^ , X 
X {i?i,.3+i(</.-(xi)|<^-(x2),£f,r)G.(x2,y) - 

- iii,.3+i(<^-(x2) I r{x{),Cf--^)Ga{x2,y)] C^{x2Zf-U{x,Zf-\ 



Finally, using (p8|) , the computation of /Fo- yields 
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FIG. 7 The partial results for a > 1 



JF^ = y dgXidgX2 Gaix,Xi)[(l){xi)P ^,(j){x2T '^]^_^Gr{x2,y) 
A/2 

= / dgXidgX2 Ga{x,xi) ^ ^ (/)(j;i)^j(/)(j;2)i2 X 
= / dgXidgX2 Ga{x,xi) ^ (/)(a:i)^^(/)(a:2)i2 x 

X {iil,.3+l(r(xi) I r{x2),£f-^) - 

0-3! 

- i?i,.3+i(r(:r2) I r{xi),Cf^?)} Hx2)lf-'<P{xi)lf-'Gr{x2,y). 

We have subsumed the partial results graphically in figure [71 where the encircled double arrows 
depict G^-chains. It is straightforward to check the cancellations — 11^ — Ilia + I^ij = 0. This 
closes the proof of the "out" -field CCR. 



Non-quasifree representation One of the main claims of this article is that, on non-stationary space- 
times, the outgoing field is in general in the GNS-representation of a state which is not quasifree. 
The reason for this is the lack of both spectral conditions and energy-momentum conservation, 
which would assure the vanishing of higher order truncated Wightman functions of the "out" -field 
in the stationary case. 

To see this explicitely, let us consider as an example for a non-stationary spacetime M*^ with 
metric (7(e) := (1 + eh)r], where rj is the Minkowski metric and h a C°°-function on M of compact 
support. One then has dgX = -y/|g[(ix = (1 + e/i)'^/^dx and such a spacetime is non-stationary since 
the metric depends on "time". 

By means of the methods described in section IIIII (see Hac07f | for a detailed calculation) , one 
computes the truncated 4-point function of the "out"-field to first order in —A in cp'^ theory on this 
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spacetime as 



r ^ 

{n,r'''{xi)r''\x2)r"\xs)r-\x,)n)J = mm J dyllD-^^{xi,y){l + eh{y))"2, (29) 

mi/^ 1=1 



where the D~ bear the subscript (7(e) to emphasize that they depend on the metric via the Klein- 
Gordon-equation. To calculate ([29l) up to first order in e, one needs to expand (1 + e/i)'^/^ and D~ 
in e, where D~{x, y) = {VL, can be expanded in e by expanding each of the two fields 

(j)^^ separately via the free Klein-Gordon equation. Denoting the expansion of (/>™ up to first order 
in e as 



we obtain the expansion of D~ to first order in e as 

D-{x,y) = {^J^{y)<t}^{x)n) + e{{^lJ^{y)<l}^{x)^) + {^J^{yW^{x)n))+0{^) 
= ■■ DQo{x,y) + €(^DQ^^{x,y) + D^Q{x,y)) +0{e^). 

To compute the single terms in the expansion of 0™, one first evaluates the wave operator as 

a = -\g\-hbg'%\hc 



d 



-hrf'dbdc - ( 2 + 1 ) idbh) d. 



□0 + eDi 



Inserting this into the Klein-Gordon equation (with minimal coupling) for 0™, one obtains to zeroth 
order in e 

(□0 + m2) < = 0, 
i.e., the Klein-Gordon equation on flat spacetime, and to first order in e 

□iC + (Do + rr?) = 



(□0 + m2) cf>^ 



d 



^0 



(30) 



Gr 



r,0 



d 



^0- 



with Grfi denoting the retarded Green's function on Minkowski spacetime. 
The well-known expression for D~ on fiat spacetime is 



Don{x,y) 



dk 
2^ 



(31) 



from which it is explicitly seen that the Fourier transform of DQQ{x,y) w.r.t. x, respectively 
X — y, has support in the negative mass shell and its Fourier transform w.r.t. y has support in the 
positive mass shell, i.e., Dqq{x, y) fulfils the spectral condition. From (f30l) it follows that D^q{x, y) 
{Dq-^{x, y)) can be obtained by application of the operator Grfi [(f + l) Vj""^ (dah) db — hm"^] to the 
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first (second) argument of Z?(^g(x,y). Fourier transforming DQ^{x,y) on Minkowski spacetime 
w.r.t. y, one thus gets the Fourier transform of Grfi multiphed by the Fourier transform of (the 
derivative of) h convoluted with the Fourier transform of (the derivative) of Dqq. Since the latter 
convolution smears up the mass shell spectrum of (the derivative of) Dqq and Gr,o is known to have 
off-shell spectrum, the Fourier transform of DQ^{x,y) w.r.t. to y clearly has off-shell support. In 
contrast, the Fourier transform of D^Q(x,y) w.r.t. to y remains on-shell like the Fourier transform 
of DQ Q{x,y). 

With these considerations in mind, we can continue examining (I29p . An expansion up to first 
order in e yields a zeroth order term and three different first order terms, viz., 



Im / dy []D^-^)(xz,y)(l + e%))2 
1=1 

4 ( d r ^ 



1=1 



+ 



j dyDQ^^{xj,y)WDQ^Q{xuy) + Y j dy D^^Q{xj,y)J\_DQQ{xuy) 

j Dd 't^J i rod 



+ 0(e2). 



(32) 



Upon Fourier transforming w.r.t. y, it is easily seen that the zeroth order term vanishes due 
to energy-momentum conservation and the spectral support properties of Dqq. Similarly, owing 
to the spectral support properties of D^q, the last first order term in ()32p also vanishes due to 
energy-momentum conservation. However, the remaining two first order terms are in general non- 
vanishing: regarding 



/ dy Y[DQ^Q{xi,y)h{y), 
L 1=1 



we can see that it does not vanish in general, despite the spectral properties of Dqq, as energy- 
momentum conservation is violated due to h not having "(5-support" in momentum space. In 
contrast to this. 



dyZ)i o(xj,y) J|L>o o(xi,y) 



is non- vanishing because of the spectral properties of D^q, even if energy- momentum conserva- 
tion holds. It might be possible to fine-tune the situation in such a way that the two above- 



mentioned non-trivial contributions to ($7,0 (xi)0 (x2)0' 



{x4)VL){ due to alleviated 



spectral properties on the one hand and abolished energy-momentum conservation on the other 
hand cancel each other, but in general this is certainly not the case. 

To close this section, we remark that from the discussion of the above example it follows that 
the metric and hence the curvature of spacetime, must show characteristic changes on a time 
scale t ^ l/(4m) to allow the violation of energy-momentum conservation (and presumably also 
the deviation from the spectrum condition) to be big enough to assure a non-quasifree "out" -state, 
e.g.,t<7x W-^s for a pion with m ~ 140MeV. This time scale is sufficiently shorter than the 
period of nucleosynthesis at about 1 to 10^ seconds after the Big Bang, such that these findings 
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are not in contradiction with well established physical facts. It is however significantly longer than 
the time of, e.g., electro-weak symmetry breaking, which happened at ~ 10~^^s such that it seems 
highly reasonable that full curved spacetime QFT calculations are required to model the physics 
of the very early universe. 

VII. UNITARY TRANSFORMATIONS BETWEEN CCR REPRESENTATIONS 

In the previous section we have seen how non-quasifree states (for free fields) naturally appear 
in scattering theory on non-stationary curved spacetimes. If one is interested in a scattering picture 
in terms of particles, one thus needs a way to calculate the particle content of a non-quasifree state, 
i.e., a unitary transformation relating the GNS-representations of a non-quasifree and a quasifree 
state, but of course the ability to calculate such a transformation is also interesting in its own 
right. Hence, given a quasifree representation of the CCR with Fock-space J-, in this section, we 
calculate the particle content of non-quasifree representations that are unitarily equivalent to the 
given quasifree one. To arrive at such a result, we shall use the language of Wightman functionals 
and ★-product calculus, c/., appendix |A] for a short introduction and notational conventions. 

To start, let (p{x) be the operator valued distribution fulfilling the Klein-Gordon equation and 
the CCR that has been obtained via the GNS-construction from some Wightman functional (not 
necessarily corresponding to a quasifree state) W_' with GNS-Hilbert space TicNS and GNS-vacuum 
state ^'o G Wgns- Furthermore, let S^{x) be the operator valued distribution obtained from a 
quasifree Wightman functional W_ via the Fock construction given in section |TI] with ^1 £ the 
Fock vacuum. A relevant application is of course ip = (p"^^ and ^ = </>™. We assume unitary 
equivalence of both CCR representations in the following technical sense: let U : J-gns J' he 
a unitary transformation such that Uip{f)U* = ^(/) V/ € ^{M) and let ^' := C/^'o £ ^ such 
that ^' is in a dense core of some closure of the Fock creation and annihilation operators a{tjj) and 
0'\x)i "0 £ "^^jX £ . It is furthermore assumed that, for any vector T in this core, the vectors 
a^{il)i) ■ ■ ■ a^{Tpn)^ are jointly continuous in the tpi w.r.t. the (H^)^^ and the topologies, where 
a" stands for either a or a^. 

To determine U, it is enough to calculate ^' since U (p{fi) • • • 93(/n)^o = Cifi) ■ ■ ■ Cifn)^, fi G 
'D(M), can be calculated using ([7]), once the n-particle components of ^ are know. Being an 
element of J^, ^ can be parameterized as 

oo „ 

^ = ^ / dgXi---dgXnfn{xi,...,Xn)C{xi)---C{Xn), 

where the complex functions are symmetric under permutation of arguments, purely positive fre- 
quency, i.e., S®^ fn = 0, and fulfil a normalization condition, viz., ||^ I'jp = Y^^=q ll'5'^"/nll(-^+)i8i" 

= 1. Furthermore, the fn are taken from some function space s.t. /„ i— > S^'^fn G {Ti'^)®'^ is 
onto, where denotes the symmetric tensor product. Given W_ and W , computing U is hence 
equivalent to determining /„, or rather S®^ fn, since the mapping S+ is not one-to-one and only 
the solution part of is "visible" in 
Let / = (/o, /i, /2, . . .), then obviously 

oo 

W: =3rDiW= Df*Df^ W, (33) 

n,j=0 

where the convergence of the infinite sums on the right hand side follows from the assumption that 
\I' is in a core for the closed creation and annihilation operators. The operators Df = Yl^=o Df„ act 
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by inserting /„ into the first /last n arguments of a Wightman function, c/., appendix[Xl Application 
of the relation (jA2h derived in that appendix and ★-multiplication with W*~ = exp^ 

{-W.^) yields 



exp 



oo „ 



E 



I6P{{l,...,n+j} 
/ = {/l,...,/fe},fc>l 



X fn{xi,---,Xn)fj{Xn+l,---,Xn+j). (34) 

Both W_ and W_ induce CCR representations. By [CoTaOSl . Lemma 5.2], this is equivalent to 
ImWg"'" = \D S'^d W^{xi^ . . . , Xn) being symmetric under permutation of arguments ~ and hence 
real by the Hermiticity of - for n > 3. This automatically holds for the quasifree state W_ 
since W2 = and W^{xi, . . . ,rc„) = for n > 3. Hence, the left hand side of is real and 
symmetric. We note that 







for > 3 



Dii = {D+{xj^ , XjJ, 0, . . .) for k = {ji , ia}, ii < 32 



Di.Wj = {0,D+{xj 



5 y ? '-'5 



0, 



for /; = {j} 



Di, = (0, Z)+( • , X,), 0, . . .) for Ii = {j} 



(35) 



As a result, only the partitions that consist of sets with one or two elements only contribute in 
Given such a partition / = {/i, . . . , let S := \Ji.\j^\=iIi and / G V'{{1, . . . , n + j} \ 5}) 
the remainder, which is a pairing partition. Employing this notation, we can can compute 



E 



♦f=i D,, w2 



E 



JeP({l,...,n+j} 

J={Jl,...,i-fc},fc>l 



5C{l,...,n+j} 
/G-P'({l,...,n+m}\5) 

/={/l,...,/(„+j_|S|)/2} 

-fi={«; Ji},«i<i! 



(n+j-|5|)/2 ^ 



(36) 



Clearly, (*jg5 D^^, ]y^)^ = if s / |5| and for s = S = {ji, . . . , jj, ji < ja < • • • < jg < n < 
ig+i < • • • < js, 

>iesBa;,]£^) {yi,...,ys 



7re5s ;=1 «=i3+l 



f^ip 2/71(0)' 



(37) 



where 5„ denotes the permutations of {1, ... Inserting ([36|) and ([37|) into ([35]) yields 



exp, - W:^) (yi,...,y.) 



E jZ j dgXi-'-dgXr E E E X (38) 

,_?;vcn "^^JW {ii,...,j4C{l,...,r} l67"({l,...,n+m}\5) 

jl<-<jq<n< i={h,..;I{r-s)/2} 

<3q+i<-<is Ii={ii,ki},ii<ki 
{r-s)/2 q s 

X n ^^^^iv^h)WD^{^ivy-n(l)) n ^'(^ivy-n(l)) X 

«=1 Z=l l=q+l 

X /n ("^l ) • • • ) X-n) fr—n (^^n+l j • • • ) 2^r) • 
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We note that fj^jdgxD^{x,y)f{x) = if / is positive/negative frequency, c/., ([7]). Furthermore, 
by our assumptions, /* is purely negative frequency and fn-r purely positive frequency. One can 
thus replace all propagators in ()38p by the real symmetric function D = + D~ since the 
integral over the added propagator with /* or fn-r always vanishes. 

Having done so, we can commute the sums over n and over S on the right hand side of (j38p . 
such that the integral contains the function 

r 

Zj. (xx , . . . , ) . — ^ ^ fn (^1 ) • • • ) Xji) fn~r (^n+1 j • • • ; Xj.) . (^9) 
n=0 

Next, we would like to symmetrise this expression, viz., 

Zr{xi, . . . , Xr) := (r!) ^ ^ Zf. (Xtt^i) ) • • • ) ^7r(r) )) 

TT&Sr 

a procedure which also makes Zr a real function, and to replace Zr by Zr- To see that this is well- 
defined, let 1 < j < r. We have to show that Zr is integrated w.r.t. a function which is symmetric 
in Xj and Xj+i. Given one term in the combinatorial sum, suppose that j,j + 1 £ S. Then, 
symmetry follows from summation over 5^. Next, suppose that either j or j + 1 is a member of a 
pairing and the other index is in S. Then, there exists another contribution to the combinatorial 
sum where j and j + 1 are exchanged showing symmetry for this case. Let finally j and j + 1 both 
be members a pairing. If the pairings are different, the argument just given applies. If this is one 
and the same pairing, then symmetry follows from the symmetry of D. 

Taking into account that the sum over Sg yields a factor si, the sum over S a factor and 
the sum over pairings a factor 2'^~'"(r — s)!/((r — s)/2)!, one obtains a combinatorial factor Cg^r by 
multiplication of these contributions. These considerations finally lead to 

oo (r-s)/2 
eWi,{W^^-W.^j iyi,---,ys) = ^ Cs,r / dgXi---dgXr Y{ D{x2l-l,X2l) X (40) 

^ r — s ^ ] 1 

Mr t — 1 

_ _ . ___ 

r 

X Yi D{xi,yi-r+s)Zr{xi,...,Xr). 

l = T~s + l 

To fulfil our task to compute U , we need to solve this system of equations for the solution part 
of Zr, i.e., for S'^'^ Zr = D'^'^Zr- To this end, to obtain a better understanding of the structure 
of (|40|) ■ we introduce some additional graphical notation: we denote the s-point function of the 
functional on the left hand side by a white circle with s legs and the function Zr by a shaded circle 
with r amputated legs. The integrations with the propagators D then either add free legs that 
carry two arrows with opposite direction or a line of that type that goes back into the shaded circle. 
S^^Zr thus corresponds to a shaded circle with r legs with double arrows of opposite direction. 
This way, one obtains two decoupled systems, one for s even and one for s odd, c/., figure IVTll for 
the even system, which makes the upper triangular structure visible. In the following, we focus 
on solving the even system, the odd system can be solved alike. In (/>P-theories with p even, the 
odd system is identically zero on the left hand side and hence gives S^'^Zr = for odd r. We note 
that the empty circles are solutions of the Klein-Gordon equation in each of their legs. By the 
demand of continuity of the creation/annihilation operators in ip and x when repeatedly applied 
to ^, Riesz' lemma implies that the empty circle with s legs is in TC^^. Let {hj}j^fq be an ONS 
in Ti.. Taking the scalar product with hj in the first two legs and then summing over j on the 
right hand side induces an opposite double arrow line that goes back into the shaded circle, since 
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1 = o 



= co,o 



+ C0,2 



+ C0,4 




+ C0,6 



+ 



0=: = 



ciam^ + C2,4 . 




+ C2, 




+ ... 




C4,4 




+ C4. 




FIG. 8 The triangular system of equations for the functions S'^'^Zr 



^ = do,o Q + do,2 Qd + do,4 + do,6 




+ 



d2,2QC^ + d2,4(^^^ + d2,( 






+ d4. 




FIG. 9 Solution to the triangular system 



{Df, Dh) = D{f, h) for /, h real. On the left hand side, we denote this contraction operation by a 
arrow-less line going back into the white circle. 

The unique solution of the even system may hence be written down in graphical form as in 
figure [9l The solution exists by assumption of unitary equivalence, thus, all infinite sums involved 
in the inverse system converge, which follows from lim„^oon(n)^ = ^ in J^, where n(n) projects 
on states with at most n particles, and the fact that for a state with at most n particles we have a 



are defined as the entries of the inverse 



finite system of equations. The constants ds^r 
of the upper triangular matrix C := {cs^r)s,re2N- 

Let us recall that the functions Zr have been a convenient intermediate tool, and that our 
ultimate aim is to determine the solution part of the (purely positive frequency) functions fn- 
Hence, it remains to reconstruct S®"^ fn = iS+^^/n from the functions S®"^ Zr- To achieve this, let 
us first suppose that = |/oP 7^ 0. As the state ^' is only determined up to a phase, one may 
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assume /o > 0. Then, by ([39]) . 

Sfzrixi, ...,Xr) = foSf''fr{xi, . . . , Xr) , for r G N. 

If S^"^ Zr = for r < ro and tq is maximal, rg must be even. It follows that S^'"'fn = for n < 

ro/2. Hence, there exist yi, . . . , yro/2 ^ M such that 5_ ^ O 5^ ^ z^^^yi, y.^^^/^, yi, . . . , ^^0/2) = 

\S®~ fro/2{yi,- ■ • ,yro/2)P > 0. We may fix the phase such that 5®~/ro/2(yi, • • • ,yro/2) > and 
we obtain the solution part of /„, n > ro/2, via 

5? 2 . . ,y^g/2,2;i>- • • = 5+ ^ /^^^/alyi, • • • ,yro/2)5+"/n(a:i,. . . ,x„), 

which closes the looked-for computation of U . 

To obtain a complete description of the scattering process on non-stationary spacetimes in 
terms of particles, we have to assume that the spacetime under consideration is asymptotically 
flat^ in the remote future and past, such that unique preferred quasifree states are available both 
for the incoming and the outgoing field. Then, there are associated Fock spaces, say J-\n = ^ and 
Tout-, and a combination of the scattering theory described in the previous sections and the results 
obtained in this section gives the n-particle amplitudes fn of the scattered incoming quasifree state 
in the particle picture of the remote future. If one wants to determine particle production from 
an incoming multi-particle state, one can apply suitably smeared incoming fields (/>™(2;) to the 
incoming vacuum, then calculate the outgoing representation of the CCR, and then conclude as 
above. 



VIII. CONCLUSIONS AND OUTLOOK 

In the this work, we have seen that non-quasifree states for free fields appear naturally in 
sc attering theory on non-stationary curved spacetimes. This result is well in line with recent works 



)ry 

HoRu02l . ISan09l | which show that a certain class of non-quasifree states, namely, the ones for which 



the truncated 2-point function is a distribution with the singularity structure of the Minkowski 
vacuum state and the other truncated n-point functions are smooth, is the natural class of states 
in perturbative quantum field theory on curved spacetimes. In the light of this, it seems somewhat 
unnatural and unnecessary to restrict oneself to quasifre e states, althoug h some important technical 
results are only available for quasifree states, see, e.g., [LiiE,o9nl . lvCT92l ]. 



Therefore, and also because there are situations where one is interested in the particle inter- 
pretation of non-quasifree states, we develop a method to calculate, provided it exists, a unitary 
transformation relating a non-quasifree state to a quasifree one. Heuristically, the form of the 
result could have been anticipated: as we assume unitary equivalence, the GNS-vacuum associated 
to the non-quasifree state corresponds to a state in the Fock space related to the quasifree state 
under consideration, and the task is to compute the n-particle components /„ of this state. Since 
we assume both states to fulfil the same commutation relations, they only differ in the real and 
symmetric part of their 2-point function and the higher order truncated n-point functions, which 
are real and symmetric in the non-quasifree case and vanishing in the quasifree case. It is thus not 
surprising that our result ()40p relates the truncated n-point functions of the non-quasifree state 



^ Here, asymptotically flat is meant in a rather loose sense, i.e., we assume that both in the remote future and in the 
remote past of {M,g), there is an open, non-empty, globally hyperbolic subset of (M,g) which contains a Cauchy 
surface of (M, g) and is isometric to such a subset of Minkowski spacetime. In this setting, it is straightforward 
to define preferred states as the pull-backs of the ones in Minkowski space. However, even within t he more strict 
definition of asymptotically flat spacetimes, one can obtain preferred states, as devised in [DMP06| . 
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to an expression in the symmetric part of the 2-point function of the quasifree state smeared with 
a real and symmetrised version of the /„. The non-trivial part of our result are, however, the 
combinatorics involved, and we have managed to tame them by encoding them conveniently into 
★-calculus on the dual of the Borchers-Uhlmann algebra. 

The method of computing a unitary transformation relating the GNS-representations of non- 
quasifree and quasifree states developed in this work is well-suited for general treatments of the 
topic, but not for explici t numeri cal calculations. A different method to compute such a transforma- 
tion, which is based on [CaG169[], works for finite-dimensional systems, i.e., "mode-by-mo de", ari d 
is therefore better suited for numerical computations, has been developed and applied in Hac07 |. 
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APPENDIX A: ★ -calculus 

In this appendix, we present some useful applications of ^-calculus and the related notation we 



have us ed in the main body of this work, ^-calculus goes bac k to Bor chers Bor75l | and Ruelle 



Rue69l |: for ★-products in algebraic quantum field theory, c/., [bDF09| | and the references cited 



therein, f or ★-calculus for quantum fields in the context of Hopf algebras see, e.g., Mestre and Oeckl 
MeOeOd ]. 



Let P be the Borchers-Uhlmann algebra of the Hermitian scalar field with multiplication (8) 
and unit 1 := (1,0,...), i.e., the non-commutative, unital, involutive, topological, free tensor 
algebra over the space of complex valued test functions T){M). If / G P is a monomial, / = 
(0, . . . , 0, /n, 0, . . .) then we identify / with /„ G 2?(M"). We also note that the involution * 

acts via /*(xi, . . . , x^) = fn{xn, . . . , xi). For fj £ V{M) and A'' C N finite, we define (dj^wfj 
{'^$fj •= 1) such that the tensor product preserves the natural order of N. A co-commutative 
co-product A : P — > P Cg) P can then be defined by 

A{0neNfn) = ^ (8>ne5/n) ® {®n&N\sfn} , (Al) 
SCN 

linearity and continuity. Note that, in (|Aip . the tensor products in the parantheses are multiplica- 
tion in P, whereas the tensor product between the paratheses is the one in P P. Furthermore, 
the projection e((/o, /i, . . .)) = /o defines a co-unit. 

Let P' = {{Wo,Wi,...) : Wo G C, W„ E P'(M^"),n > 1} be the topological dual of P, i.e., the 
space of Wightman functionals. In many applications, a Wightman functional will be given by the 
sequence of n-point (truncated) VEVs, viz., wji"\xi, . . . , x„) := {Q, ip{xi) • • • ip{xn)^)^'^^ for some 
operator valued distribution (/?. The co-product on P naturally leads to a product, ★, which can 
be defined as 

Wi.V := {W0V)oA, 
making a commutative algebra with unit 1 = e. 



From ()Aip we get 

iW*y)\N\{'S)neNfn) = (®nGs/n) V|Ar|_|s| ((g)„e7v\5/„ 



SCN 
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which shows the comcidence of * with Borchers' s-product. 

It is easy to see that {W*'^)m = for n > m and W G V^^ : = {Z S P' | Fq = 0}. Hence, 
arbitrary *-power series converge on P']^. In particular, the ^-exponential exp^ : P'^ ^ 1 © P'^ 
and logarithm log^ : 1 © P'^ — > 'D^i are well defined through their power series exp^(l^) : = 
'^'^=oW*"' /n\ and log(VF) := —Yl'^=ii^ ~ W)*^/n. Furthermore, expj, and log^ are inverses of 
one another and the usual relations hold, i.e., exp^fM^ + V) = exp^ W * exp^ V, exp^(O) = 1 and 
log^fiy -k V) = log^H^ + log^y. For log^ W, we also use the notation . In fact, if W_ denotes 
the collection of Wigh tman functions, (W'^)n is the truncated n-point function defined in equation 
© above, c/., [Bor75l |. 



Let / G 'D{M), and m {f),m (/) : V ^ V he the left and right multiplication with / in 

P. We then define left and right derivatives Dj and Dj on P' via Df W_ := W o m (/) and 

Df W_ := Wo m (/). In fact, it is easily verified that a Leibniz rule holds for these derivatives 
(which of course motivates this nomenclature), viz., 

Df (WkV) = (Df W)irV + Wk{Df V) 

where Df stands either for Df or for Df. For a formal power series h{t), this implies that Df 

h-k{W_) = (h')^,(W)-k Df W_, where denotes the formal ★-power series one obtains from the formal 
power series / by replacing the normal product with the ^-product and f'{t) := f^^\t) := df/dt. 

For notational convenience, we shall sometimes write Df=. J^.jdgxf(x) Dx- 

Higher order derivative operators Df W_ '■= W o m (/), / G P, can also be written as Df= 
E~=o Df„, and we write 

Df„ = : j dgXi---dgXnfixi,...Xn)Dxi---Dx„- 

M" 

By induction, it is easy to see that, for and fj symmetric, the following chain rule holds 

DfnDf^ K{W) = ^{h^%{W) ★ / dgXl ■ ■ ■ dgXn+J *'=1 ^ ^ 

k = l /6P({l,...,n+,}) 



X fn(,Xl, . . . , Xn'jfj {Xn+1 j • • • ) Xn+j ) , 



where V{N) is the set of partitions of and Dii ■ =Dxj^ ••• Dxj^Dxs^ "■■ Dxsr ~ 

{jl, . . . ,jq,Si, . . . , Sr}, I < jl < j2 < ■ ■ ■ < jq < n < Si < S2 < . . . < Sr < n + j. 
Fnr h(f^ = pvTif^V W"^ P T)', and W = pyn (W^\. wp finallv nhta.in 



For h{t) = exp(t), W^ G V^^, and W = exp^(W:^^), we finally obtain 

B/„B/, W = W* / dgXi ■ ■ ■ dgXn+j Yl ^ ^ 

, , N -fST'fll,. ..,n+j} 

Mx("+^) /={/!. ..,/^},fe>i 

X fn{xi, ■ ■ ■ , Xn) fj{Xn+l, • • • ,Xn+j), (^2) 

which closes this appendix. 



APPENDIX B: Identities involving retarded products 



In this appendix, we would like to prove some identities we have used both implicitely and 
explicitely in the proof of the CCR of the outgoing field in the main body of this work. 
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First, we need to show that the cr-th order of the j-th power of the interacting field is, hke in 
the case of the interacting field itself, the result of the cr-fold retarded action of £int on the j-th 
power of the "in" -field. For this purpose, it is convenient to introduce a shorthand notation for 
the multiple application of the formal differential operator V^x), viz., 

'D'^A{x) := J dgXn,iT>'^{Xn,i)A{x)lx„^i^x ■= J • • 



where X^^i denotes the causally ordered /c — / + l-tuple {xk,Xk-i, • • • , xi)lx,.^-yxi and A{x) is an ar- 
bitrary expression in terms of "in"-fields. Later, we will also need tuples of spacetime points devoid 

of mutual causal relations. To this end, we will denote by the k — l + 1-tuple (x^, . . . , x/). 
With this notation, we have to show that 

J]o-j=(T« = l 

On can achieve this by an induction similar to the one employed to show that (j){x)cj defined via 
the Yang-Feldman equation is equal to P°"</)™(x). The case o" = is clear. Let us assume the result 
holds for fj and show how its validity for o" + 1 follows. We can compute 

v''+^(l}'\xy = vv''(i}'\xy 



lW{x)lx^^^yX^^,^^x 



: J dgX^+iV{x^+i) n ( / dgX^.A^^'^'i^a.. 

dgX„^^^V">'+\x„+i X X.„i)0'°(x)l,^+,^x.,,,^,x (Bl) 





In comparison, we know that (/)(x)^^_-^ equals 

dgX.,,iP-«(X^,,i)<A-(x)lx.^,,^,j , (B2) 

notably, the appearing j integrands in every summand of X^j^o- =(t+i ^'^^ mutually independent, 
while the corresponding integrands in ()Bip depend on one another, namely, only one of them 
depends on x^^j^i explicitely, but all others are constrained such that their integration domains 
are causally later than x^+i . If we can show that all combinatorially possible constraints appear 
exactly once, then their sum yields the independent integrands of (jB2|). To achieve this, let us 
pick an arbitrary but fixed summand of YliY,<Ti=cT+i^ ^^y-> ('^i' ■ ■ ■ ^^j) — (^ii ■ ■ ■ ,nj). Apart from 
integration domain restrictions, the j integrands of this summand are the same as the ones of all 
summands of Ylj2'^i='^ Sfc=i ™ (IBlh with (cti, . . . , 0"^ + 1, . . . , (Tj) = (ni, . . . ,nj). If we assume that 
m entries of (ni, . . . ,nj) are non-zero, then there are exactly m summands of the latter type due to 
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the exhaustion of the sum and the Leibniz rule for 'D{xa+i)- The integrands of these m 

summands are such that it is always a different of the j integrands that explicitely depends on x^+i- 
Therefore, all possibilities of the event "One of the j integrands involves an integration variable 
which is causally earlier than the integration domains of all other j—1 integrands." appear, and their 
sum yields exactly the summand of ()B2p under consideration; this implies V^^^ cp^'^ {xy = (j)(xy^^-^. 
In the following, we can use this result to prove the looked-for recursive relations 



d,x, Gr{x,x,) Y.^i^^y'^. {rixi) I r{y),ci?) €f~Kxi) m 



M 



p-2 

E 

J2(Ti=n-lj=0 
p-2 



To this effect, let us define another formal differential operator 6{x) via 

5{y)A{x) :=D{.,y)^0^^ := ^ [r{y), A{x)] , 

where A{x) is again an arbitrary expression in terms of incoming fields and the empty argument 
(• ) denotes that the argument of the "in"-field the differential operators currently acts upon in a 
Leibniz rule-summand has to be inserted in that slot. 

We can now proceed to prove (IBSh by induction. In case n = 1, one can straightforwardly 
compute that both expressions equal 

(p-l)i j dgXi Gr{x,xi)Gr{xi,y)(l){xi)P~'^. 

M 

Assuming that the wished- for equality holds for n, we can analyse the case n + 1, viz., 



n+l n+1 



n + l ^^ = 1 

„ n+l n 

n + l « = 1 k = l 

+ i6{y) J dgXn+1,1 P"+H^n+i,i)0'"(a;)lj,^x„+i,ibx 

/n+l , -Nn+l ^/ 

n dgy^ V{xn+iY~^R+n+i i^ix) \ (y) , £f,^ (x£ )) 1 . 

n+l -'- ' 

+ i5{y) j dgX^+i,! G^(x,xi)P"(X„+i,2)(/>^°( XlY '^yyXn+\^2>ixiysi 
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The remaining steps are cumbersome but elementary. To write them down at this point, we 
would have to introduce even more abbreviating notation which is necessary to avoid loosing 
the overview, we thus prefer to rather sketch them briefly: the next step would be to insert the 
induction hypothesis in the first term of the last line of the above equation and to rewrite the 
resulting retarded product of lower order by means of the formal differential operators D and 
5. Employing the Leibniz rule for T> and the combinatorial considerations regarding different 
possibilities to distribute causal relations among the factors of a product we have already used 
at the beginning of this appendix, one can show that the two terms in the last line of the above 
equation indeed add up to 

which proves the first of the two looked-for recursion relations for Ri.n+i (i?5>"(x) | </>™(y), . 

To prove the second one ()B4p . we can again perform an induction by order. The case n = 1 can 
be validated straighforwardly and for n + 1 a computation yields like in the preceding proof 

„ n+1 n 

= i J ]ldgX,V{Xn+l)J2{^''~'(^nMiny)VHXk,,)r{x)^^^^ + 

+ i5{y) J dgXn+1,1 P"+i(X„+i,i)0*'^(x)lj,^x„+i,ib- 

Employing the by now multiply used combinatorial considerations, one can show that the first 
summand in the last line equals 

up to some missing combinatorial possibilities. To show that the second summand in the last line 
accounts for exactly these possibilities, one has to be able to "extract" a Gr{- ,y) from it. This 
can be achieved by the Leibniz rule for 5. To wit, every V in S{y)'D^~^^{Xn+i,i)(j)^'^{x) involves a 
commutator with £int = {4>™Y /p, such that the application of 5{y) to V^~^'^{Xn+i^i)(t)™{x) yields a 
sum over terms without 5{y) but with one T> replaced by a commutator with (5(y)>Cint, which equals 
Gj.{ ■ , • due to the enforced causal relations of the arguments. Further elementary steps 

then lead to the wished-for conclusion. 
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